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Abstract 



We have developed a three-dimensional numerical model and applied it to simulate plasma 
flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic 
field. The model is based on the assumption that the plasma dynamics are determined by 
the slow mean flow, which forms a backdrop for the rapid propagation of MHD waves. The 
equations describing the slow motion of matter were obtained by averaging over rapidly 
propagating pulsations. The numerical model includes the diffusion of magnetic field by 
current dissipation in turbulent vortices, magnetic buoyancy, and wave MHD turbulence. 
A modified three-dimensional, parallel, numerical code was used to simulate the flow 
structure in close binary systems with various accretor magnetic fields, from 10 5 to 10 8 G. 
The conditions for the formation of the accretion disk and the criteria distinguishing the 
two types of flow corresponding to intermediate polars and polars are discussed. 

1 Introduction 

Among close-binary systems in which magnetic fields exert a substantial influence on the flow 
of matter in the system, the most prominent are polars and intermediate polars [1] . These are 
semi-detached binary systems consisting of a low-mass, late-type star (the donor) and a white 
dwarf (the accretor). During the process of mass transfer, matter flows from the donor onto 
the white dwarf through the inner Lagrange point. 

In polars (AM Her- type systems), the white dwarf has a strong intrinsic magnetic field (B a ~ 
10 7 — 10 8 G at the surface). These systems are characterized by relatively short orbital periods 
(1-5 hr),and the rotation of the components is synchronized [2]. Analysis of observational data 
shows that accretion disks do not form in polars; instead, the matter flowing from the donor 
forms a collimated flow, which falls onto a region around one of the magnetic poles of the 
accretor along the star's magnetic-field lines [U [3]. 

The white dwarfs in intermediate polars possess relatively weak magnetic fields (B a « 
10 4 — 10 6 G at the surface). These systems occupy an intermediate position between polars and 
non-magnetic cataclysmic variables. The rotational periods of the accretors in these systems 
are appreciably shorter than the orbital periods, by factors of one, two, or even three orders 
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of magnitude [2]. The asynchronicity of the rotation of the accretor in intermediate polars 
can be explained as an effect of the interaction of the magnetic field of the white dwarf with 
the disk material near the boundary of the magnetosphere. As a result of this interaction, 
an equilibrium rotational regime is established in which the corotation radius is equal to the 
magnetosphere radius [H H] . 

Our previous studies [5j El [7] have used three-dimensional numerical simulations to in- 
vestigate the flow structures in close-binary systems (using SS Cyg as an example), allowing 
for relatively weak, dipolar accretor magnetic fields (B a = 10 5 G at the surface). Model- 
ing with high values of B a is more difficult due to the fact that increasing the field strength 
leads to rigid constraints on the time step required for the difference scheme to be stable (the 
Courant-Friedrichs-Levy condition). However, a simple analysis shows that this problem has 
a fundamental, rather than technical, nature. Indeed, the Alfven velocity, estimated by the 
expression 



can exceed the speed of light c in the accretion flows of polars. Therefore, in general, it is 
not correct to use a non-relativistic magnetic gas-dynamical approximation when modeling 
the structure of such a flow. However, it does not really make sense to apply relativistic 
magnetic gas-dynamics in this case, since the plasma flow itself is appreciably non-relativistic. 
Only Alfven and magneto-acoustic waves in the accretion flow will propagate with relativistic 
speeds. Over the characteristic dynamical time for the evolution of the slow plasma flow, the 
MHD waves are able to pass through the accretion flow many times (in the longitudinal and 
transverse directions). As a result, the dynamics of the plasma flow can be treated in a modified, 
non-relativistic, magnetic, gas-dynamical framework, as a mean flow that forms the backdrop 
for wave MHD turbulence. 

We propose here a numerical model for computing the flow structures in close-binary systems 
with strong accretor magnetic fields. We present results of our numerical simulations obtained 
using a modified, three-dimensional, parallel code. These computations were carried out for 
magnetic fields at the accretor surface from 10 5 to 10 8 G. 

The paper is organized as follows. Section 2 describes the main equations and numerical 
method used in the framework of the modified model. Section 3 presents a comparison of the 
flow structures in a semi-detached system similar to SS Cyg for various white-dwarf magnetic 
fields, based on the three-dimensional numerical simulation results. The Conclusion briefly 
discusses the main results of our study. 

2 Numerical method 
2.1 Main equations 

We will use a non-inertial reference frame that rotates together with the binary system with 
the angular velocity ft relative to its center of mass. We will use the Cartesian coordinates (x, 
y, z) in this frame, with the coordinate origin coincident with the center of the accretor and 
the center of the donor located a distance A from the accretor along the x axis. The z axis is 
directed along the rotational axis of the system. 

Let us consider a semi-detached binary system whose accretor rotates synchronously, so 
that the rotational period of the accretor is equal to the orbital period of the system. We will 
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take the magnetic field of the accretor to be dipolar, so that the vector magnetic field is given 
by 

In general, the direction of the magnetic moment /x may not coincide with the direction of the 
angular velocity of the binary system O. In our chosen coordinate system, the vector /j, has 
the components p x = /i sin 9 cos 0, \i y = psin^sin^, p 2 = pcos9, where p is the modulus of n, 
9 is the inclination of /x to the z axis, and is the angle between the x axis and the projection 
of \i onto the xy plane. Note that a magnetic field of the form ([2]) is potential V x B, = and, 
in the case of synchronous rotation, stationary dH*/dt = 0. 

The magnetic field of the accretor can be fairly strong in the region of the magnetosphere. 
Therefore, it is convenient to represent the total magnetic field in the plasma B as the sum of 
the accretor magnetic field B* and the magnetic field induced by currents in the plasma itself 
b: B = B* + b. Using the properties of the field B* noted above, it can partially be eliminated 
from the magnetic gas-dynamical equations [H [91 ITU] . 

We described plasma flows arising due to mass transfer in the close-binary system, taking 
into account the strong magnetic field of the compact object B*, using the system of equations 



dp 
dt 



+ V • (pv) = 0, 



dv , t— 7\ VP b x (V x b) n . 

+ v ■ V v = \ + 2 v x VI + g 

dt p Anp 



V_L 

T 



pT 



dh 

~dt 
ds 
dt 



Vx[vxb + vxB*- 7?(V x b)] , 



+ (v- V)s 



n 2 (r-A) + ^-(Vxb) 2 . 
Air 



(3) 

(4) 
(5) 
(6) 



Here, p is the density, v the velocity, P the pressure, s the entropy per unit mass of the gas, 
n = p/m p the number density, m p the proton mass, rj the magnetic viscosity, and g = — V$, 
with $ being the Roche potential. The effects of radiative heating and cooling are included 
in the entropy equation dHJ), as well as heating of matter due to current dissipation (the last 
term). The radiative heating and cooling functions T and A have complex dependences on 
the temperature T [TTJ [T2J [121 EH]- Our numerical model uses a linear approximation for these 
functions in the vicinity of the equilibrium temperature, T = 11230 K [J51 El [TB], corresponding 
to an effective temperature of the accretor of 37000 K. The term 2 (v x ft) in the equation of 
motion describes the Coriolis force. The density, entropy, and pressure are related by the 
equation of state of an ideal gas: s = c\/ln(P/p 7 ), where cy is the specific heat of the gas at 
constant volume and 7 = 5/3 is the adiabatic index. 

The numerical model includes the effects of magnetic-field diffusion [see (jSJ)]. The earlier 
analysis of [51 [6] showed that two effects dominate in accretion disks forming in close-binary 
systems '3]. The first is magnetic reconnection and current dissipation in turbulent vortices. 
The second is buoyancy of force tubes of the toroidal magnetic field generated in the disk 
due to its differential rotation. Accretion disks may not be able to form in systems with 
strong magnetic fields, in which case the flow will have the form of an accretion-column stream 
directed from the inner Lagrange point Li toward one of the magnetic poles of the accretor (a 
polar). In this case, the main effect leading to current dissipation in the plasma is probably 
wave MHD turbulence due to the propagation of Alfven and magneto-acoustic waves in the 
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accretion-column stream. The velocity of these waves will far exceed the velocity of the plasma 
itself and, in some cases, may even be relativistic. The wave magnetic viscosity (see Appendix 
B) rj w can be obtained from the expression 

l w B^ 

Tj w Oi w , yl) 

where l w is the characteristic scale for wave pulsations, which can be estimated as the scale for 
inhomogeneities of the accretor magnetic field: l w = B*/\VB )t .\. The parameter a w determines 
the efficiency of wave diffusion. The total magnetic viscosity 77 due to all effects depends on 
the magnetic field in the plasma. Therefore, on the whole, the diffusion of the magnetic field 
has a non-linear character. 

The last term in the equation of motion (j3J) describes the force exerted on the plasma by the 
accretor magnetic field, which influences the component of the plasma velocity perpendicular 
to the magnetic field lines Vj_. The basis for this term is provided in Appendix A. Note that the 
motion of plasma particles across the magnetic field is due primarily to the gravitational force 
of the compact object (gravitational drift) (see, for example, [El HH HH]). Due to the Larmor 
character of the motion of particles in the magnetic field, their mean motion is decelerated in 
the perpendicular direction. The strong external magnetic field plays the role of an effective 
fluid with which the plasma interacts. The last term in can be interpreted as a frictional 
force between the plasma and magnetic field, whose form is analogous to the friction between 
the components in a plasma consisting of several types of particles [T7j. Recall that our model 
assumes synchronous rotation of the accretor. In this case, the velocity of the magnetic field 
lines in our chosen system is zero. The characteristic time for decay of the perpendicular 
velocity is 

4irpr] w 



Bi 



(8) 



This quantity is determined by wave dissipation of the magnetic field, which is characterized 
by the diffusion coefficient r] w (J7j). A similar expression for the electro- magnetic force due 
to the strong magnetic field of an accretor was used in [201 ED E21 E31 El EU E5] , where the 
flow structure was modeled using a quasi-particle method. The quasi-particles were taken to 
be individual blobs of plasma that fall into the Roche lobe of the accretor through the inner 
Lagrange point L\. In these studies, an expression for the decelerating force of the form — v^/r 
was justified by the fact that diamagnetic effects arise in the plasma blobs during their motion 
in the external magnetic field |26] , which hinder their free motion perpendicular to the magnetic 
field lines. 



2.2 Numerical method 

Suppose we know the distribution of all quantities in the computational region at time t n . To 
obtain their values at the next time step, corresponding to time t n+l = t n + At, we split the 
equations according to the various physical processes they represent. The full algorithm 

consists of five successive steps, which we describe briefly below. 

In the first step, a subsystem of equations describing the dynamics of the plasma in its own 
magnetic field is distinguished: 

^ + V-(pv) = 0, (9) 
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|? = Vx(vxb), (11) 

- + (vV), = 0. (12) 

The form of this system coincides with the equations of ideal magnetic gas dynamics [27J. 
The system can be solved numerically using the higher-order Gudonov-type difference scheme 
described in [T6] (if B* = in all expressions there). When constructing the difference scheme, 
we applied the technique of unified variables for MHD equations [2B] , which enabled us to use an 
adaptive mesh in the numerical code. The computations presented below used a geometrically 
adaptive mesh that became more dense toward the equatorial plane and the surface of the 
accretor. This made it possible to appreciably enhance the spatial resolution of the vertical 
structure of the accretion disk and in the region of the accretor magnetosphere. We used the 
eight-wave method to clean the divergence of the plasma magnetic field b [HI • Note that 
the external magnetic field B* does not appear in this system of equations. Therefore, the 
limitations on the time step At imposed by the stability condition for the difference scheme 
(the Courant-Friedrichs-Levy condition) will not be excessively severe. 

In the second step, we take into account variations of the gas velocity due to external forces 
(the Coriolis force and gradient of the Roche potential): 

— = 2(vxft) + g. (13) 

All remaining quantities are taken to be constant in this step of the algorithm. Since the Roche 
potential is time- independent, the solution of this equation in the interval t n <t < t n+1 can be 
written 

v x = A cos(2fit) + Bsin(2Qt) + J^, (14) 

v y = Bcos{2ttt) - Asm{2Qt) - (15) 
v z = v ° z -g z (t-t n ). (16) 

The integration constants in these equations A and B are determined from the initial conditions: 

A = ( v * ~ In) cos ( 2 ^") - {< + m) sin ^")' 

B = ( v * - In) M2nn + ( v ° + 2%) cos ( 2fir )- ( 18 ) 

The initial velocities v° are equal to the values obtained at time t n+1 in the previous step of 
the algorithm from the solution of fl9l)- fTT2|) . 

The third step of the algorithm takes into account deceleration during motion of the plasma 
across magnetic field lines, as well as the generation of magnetic field due to this motion. The 
corresponding equations can be written 

dv± v I (9b _ , t-. % / \ 

-flT = -T' sF = Vx(vxxB,). (19) 
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Integration of these equations yields 

v± = v°e"^ fl , b = b° + r(l-e-^)vx(v° xB t ). (20) 

In the fourth step, we take into account the effects of magnetic-field diffusion. Recall that 
the equation describing diffusion of the magnetic field is non-linear in our model. Therefore, the 
application of explicit methods to solve this equation would lead to excessively severe limitations 
on the time step. In our approach, this equation was solved numerically using an implicit, 
locally one-dimensional method with a factorizable operator (30]. We applied regularization 
of the factorizable operator to deal with the mixed derivatives that arise due to the use of an 
adaptive mesh. The regularization procedure essentially reduces to replacing the multiplicative 
operators making up the original factorizable operator with equivalent tridiagonal operators. 
The regularization parameter is determined by the maximum-modulus eigenvalue of the metric 
tensor describing the curvilinear coordinate system. To correctly take into account non-linear 
terms in the scheme, we devised an iterative process that was applied until we obtained a 
solution with a specified accuracy. In each iteration, a system of linear, algebraic equations 
with a tridiagonal matrix arises, which is solved numerically using a scalar fitting method. The 
method used to solve the magnetic-field diffusion equation is described in more detail in [3T) . 

Finally, the fifth step takes into account the effects of radiative heating and cooling, as well 
as heating due to current dissipation. These processes are described by the right-hand side 
of ([6]). We emphasize that we used semi-analytical or implicit, absolutely stable methods in 
all other steps of the algorithm (apart from the first). Therefore, no constraints are imposed 
on the time step At apart from the Courant-Friedrichs-Levy stability condition [which arises 
during the solution of (19|)- (|T2|) ]. 

3 Computational results 
3.1 Parameters of the model 

Let us consider the formation of an accretion disk in a close-binary system whose parameters 
correspond to SS Cyg as a function of the accretor magnetic field. The donor of the SS 
Cyg system is a red dwarf with mass 0.56 M & and the accretor is a white dwarf with mass 
M a = 0.97 Mq. The orbital period of the system is P or b = 6.6 hr, and its semi-major axis is 
A = 2.O5i? [32] ■ In our computations, the magnetic field at the white-dwarf surface was varied 
from 10 5 to 10 8 G, and the orientation of the magnetic axis was fixed, 9 = 30° and <fi = 0°. 

At the inner Lagrange point L 1; we specified the gas velocity to be equal to the local sound 
speed, c s = 7.4 km/s, which corresponds to a donor temperature of 4000 K. The gas density 
at L\ is p{Ll) = 1.1 x 10~ 7 g/cm 3 , and the mass transfer rate is M = 10~ 9 Af /yr. We 
specified the following boundary conditions at the remaining boundaries of the computational 
domain: density p h = 10~ 6 p(Li), temperature T b = 11230 K, velocity v b = 0, and magnetic 
field Bb = B*. The accretor was taken to be a sphere with radius 0.0125A, at whose boundary 
we specified a free-outflow condition. The velocity of the inflowing plasma was taken to be 
parallel to B*. All the matter in a cell occupied by the accretor was taken to fall onto the white 
dwarf. We emphasize that these boundary conditions were used at each individual stage in the 
subdivided numerical algorithm described in the previous section. The initial conditions in the 
computational region were density p = 10 _6 p(-^i); velocity v = 0, temperature T = 11230 K, 



6 



1000 



100 



10 



: io;* 

10 9 

< 

m 11 


y ,.•*** 

y 






,«** 


■ 






■ ■ ■ i i 





10 q 



10" 



0.1 



0.01 



10° 
B a! G 



10' 



10 c 



Figure 1: Radius of the white-dwarf magnetosphere in the SS Cyg system as a function of the 
surface magnetic field. The individual curves correspond to various accretion rates (in units 
of M Q /yr). The solid horizontal lines show the numerical boundary of the accretor surface, 
0.0125A (lower), and the distance from the accretor center to the inner Lagrange point, 0.56A 
(upper). 



and magnetic field Bo = B*. We found our solution in the region (— 0.56A < x < 0.56A, 
-0.5QA <y< 0.56A, -0.28A < z < 0.28A) on the geometrically adaptive mesh [T6] . 

The extent to which the magnetic field influences the flow structure can be determined 
using estimates of the radius of the accretor magnetosphere. We will assume that the mag- 
netic pressure is equal to the dynamical pressure of the accreting gas at the boundary of the 
magnetosphere (see jl]): 

02 ' (21) 



in 



where Vff — ^/2GM a /r is the free-fall velocity. We can find the density of the matter from the 
expression for the accretion rate: M a = A.Tir 2 pVff. Substituting B = B a (R a /r) 3 , where R a is 
the radius of the accretor, into (l2Tj) . we can find the radius of the magnetosphere: 



( B i Rl a 



12 



\%GM a M* 



(22) 



The dependence of the magnetosphere radius r m on the magnetic field B a at the accretor surface 
is shown in Fig. [TJ The different lines correspond to different accretion rates M a in units of 
M /yr. The solid horizontal lines show the specified accretor radius, 0.0125A (lower line), and 
the distance from the center of the accretor to the inner Lagrange points, 0.56A (upper line). 
In our computations, the characteristic accretion rate was of the order of lO _lo M /yr. In this 
case, the magnetosphere radius is approximately 5 — 6R a for B a = 10 5 G and roughly 20i? a for 
B a = 10 6 G. When B a = 10 4 G, the magnetosphere radius becomes smaller than the radius of 
the "numerical star", although it remains larger than the radius of the white dwarf. 
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We present results for eight computational models here, with magnetic fields at the accretor 
surface equal to 10 5 G (model 1), 5 x 10 5 G (model 2), 10 6 G (model 3), 5 x 10 6 G (model 4), 
10 7 G (model 5), 5 x 10 7 G (model 6), and 10 8 G (model 7). These models can be divided into 
two groups. The first includes models 1, 2, and 3, with relatively weak magnetic fields, in which 
accretion disks form. These models presumably correspond to the case of intermediate polars. 
The second group contains models 4, 5, 6, and 7, with strong magnetic fields and without 
the formation of an accretion disk. These models correspond to the case of polars. The 
computations for all the models were continued until a quasi-stationary regime was reached, 
when the total mass of matter in the computational domain was roughly constant (to within 
1%). In models in which an accretion disk formed, the time for establishing this quasi-stationary 
regime was about 10-15 orbital periods. In models with strong magnetic fields and without the 
formation of an accretion disk, the transition to a quasi-stationary regime was more rapid (about 
five orbital periods). We will describe the results for these two groups of models separately. 



3.2 Weak magnetic field 
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Figure 2: Flow structure in the equatorial (left) and vertical (right) planes for model 1 (B a = 
10 5 G). The distributions of the density (contours on the left and gray scale on the right) and 
velocity (arrows) are shown. The dashed curve corresponds to the boundary of the accretor 
Roche lobe. The bold curves with arrows show either the streamlines from L\ (left) or the 
magnetic lines (right). 

Figures [2HH show the flow structures in the equatorial (xy) and vertical (xz) planes for 
model 1 (B a = 10 5 G), model 2 (B a = 5 x 10 5 G), and model 3 (B a = 10 6 G), respectively. 
The distributions of the density (contours or gray scale) and velocity (arrows) are shown. The 
values of the logarithm of the density are indicated on some of the contours (in units of p(Li)). 
The left panels of these figures also show the boundary of the accretor Roche lobe (dashed 
curves) and a streamline originating from the inner Lagrange point L\. The right panels show 
the magnetic lines. 

Figures [3H7] present the three-dimensional structure of the flow for models 1-3. These 
figures show constant surfaces of the logarithm of the density at the level lgp = —4.5 (in 
units of p(Li))and magnetic lines originating from the accretor surface. The gray scale for the 
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Figure 3: Same as Fig. [2] for model 2 (B a = 5 x 10 5 G). 
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Figure 4: Same as Fig. [2] for model 3 (B a = 10 6 G). 



magnetic field displays the magnetic field strength. The rotational axis (thin vertical line) and 
magnetic axis (bold inclined line) of the accretor are also shown. 

The flow pattern obtained for model 1 (B a = 10 5 G) is in good agreement with the results 
of our earlier computations for this case carried out for a rigorous MHD model [S]. The char- 
acteristic parameters of the accretion disk, and especially its structure, are virtually identical 
for these two sets of computations. A magnetosphere forms near the accretor surface, and the 
accretion has a column-like nature. Matter in the magnetosphere moves primarily along mag- 
netic lines and falls onto the star mainly in the region of its magnetic poles, forming accretion 
columns. Vacuum zones form in the region of the magnetic equator of the accretor. This is due 
to the fact that the magnetic field prevents matter from penetrating to these regions, since the 
magnetic lines lie predominantly along the stellar surface near the magnetic equator. We can 
note some modest differences compared to our previous computations. These are associated 
with the fact that we developed the model we consider here to simulate plasma flows in strong 
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Figure 5: Three-dimensional flow structure for model 1 (B a = 10 5 G). A constant-density 
surface (lgp = —4.5, where p is in units of p{L\j) and magnetic lines are shown. The gray scale 
for the magnetic lines displays the magnetic field strength. The thin, vertical line indicates the 
rotational axis of the accretor, and the bold, inclined line its magnetic axis. 

magnetic fields. A more refined adjustment of parameters is apparently required for the case 
of relatively weak fields (in particular of a w ) in order to obtain more complete agreement with 
the more rigorous MHD model. 

The outer radius of the accretion disk becomes substantially smaller (about 0.15 A) when 
B a = 5 x 10 5 G (model 2). The effective magnetic braking and angular-momentum transfer 
are enhanced, as is demonstrated by the behavior of the streamline shown in Fig. [3j The 
magnetosphere region has become much larger. The accretion columns near the magnetic poles 
of the accretor are more clearly visible. 

Finally, the accretion disk essentially becomes degenerate when B a = 10 6 G (model 3). 
Matter is able to undergo one to two orbital revolutions before falling onto the accretor. The 
term " spiral disk" is more suitable for describing this flow structure, since the matter velocities 
involved differ strongly from their Keplerian values. The outer radius of this accretion spiral 
disk is roughly 0.1 A. It is located virtually entirely in the region of the accretor magnetosphere. 
A substantial part of the spiral disk is comprised of the accretion-column streams. This model 
essentially corresponds to the limiting case of intermediate polars. 

Note that the size of the accretor magnetosphere region obtained in the numerical compu- 
tations is in good agreement with analytical estimates obtained using fl22l) (Fig. [I]). 

3.3 Strong magnetic fields 

The MHD flow structures for models 4 (B a = 5 x 10 6 G), 5 {B a = 10 7 G), 6 (B a = 5 x 10 7 G), 
and 7 (B a = 10 8 G) are presented in Figs. IB| l9| fTUl and fTTl respectively. 

The upper panels in these figures present the three-dimensional structure of the flow. These 
show constant surfaces of the logarithm of the density at the level lgp = — 5 (in units of p{Li)) 
and magnetic lines. The gray scale for the magnetic field displays the magnetic field strength. 
The rotational axis (thin, vertical line) and magnetic axis (bold, inclined line) of the accretor 
are also shown. The images have been turned so that features of the accretion flow near the 
accretor surface are visible. The right panels show distributions of the density (gray scale) and 
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Figure 6: Same as Fig. [5] for model 2 (B a = 5 x 10 5 G). 




velocity (arrows) in the vertical plane xz. The boundary of the accretor Roche lobe (dashed 
curve) and magnetic lines are also shown. 

Analysis of these figures indicates that the flow structures for this group of models are 
qualitatively different from those for the previous cases (models 1, 2, and 3). No accretion disk 
forms in these models, and the flow takes the form of accretion- column streams. In models 4 
and 5, the streams originating from the inner Lagrange point L\ divide into two flows when 
they reach the accretor surface, with the left and right flows falling onto the North and South 
magnetic poles of the star. The northern accretion flow is more powerful in model 4, and the 
southern in model 5. This comes about because, in our formulation of the problem, the South 
magnetic pole is closer to L\ than the North magnetic pole. Therefore, the motion of plasma 
toward the South magnetic pole is energetically more favorable in a stronger magnetic field. 
This is supported by our results for model 6, in which only one accretion flow forms, ending at 
the South magnetic pole of the star. 

The magnetic field in model 7 is so strong that it virtually completely controls the flow inside 
the accretor Roche lobe. The matter is captured by the magnetic field essentially immediately 
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Figure 8: Flow structure in model 4 (B a = 5 x 10 6 G). A const ant- density surface (lgp = —5, 
where p is in units of p(Li)) is shown in the upper panel. The remaining notation is the same 
as in Fig. [5j The lower panel shows the distributions of the density (gray scale) and velocity 
(arrows) in the vertical plane. The dashed line corresponds to the boundary of the accretor 
Roche lobe. The bold lines with arrows show magnetic lines. 



behind the Lagrange point, and is directed along the magnetic lines toward the stellar surface, 
forming a powerful southern stream and weaker northern stream. Figure [1] shows that, in this 
case, the magnetosphere becomes larger than the accretor Roche lobe, and therefore partially 
encompasses the donor envelope. This model probably corresponds to so-called magnetors [I], 
which represent the limiting case of a polar whose white dwarf has a very strong magnetic field. 



3.4 Discussion 

Thus, in our computations of the flow structures in close-binary systems whose parameters 
correspond to SS Cyg, models with magnetic fields B a < 10 6 G form an accretion disk, while 
no accretion disk forms in models with stronger fields. We can present the following simple 
arguments to explain this result. In a first approximation, the behavior of the stream of matter 
flowing from the inner Lagrange point L\ inside the accretor Roche lobe can be analyzed in a 
ballistic approximation, without including effects associated with the pressure or magnetic field 
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Figure 9: Same as Fig. [S]for model 5 (B a = 10 7 G). 

[33| 134"] . This approximation is justified by the supersonic nature of the flow in the stream. A 
trajectory analysis indicates that the stream approaches fairly closely to the accretor surface 
|35j . The minimum distance between the trajectory of an individual particle and the center of 
the star -R m i n depends on the mass ratio q = M d /M a . In the range 0.05 < q < 1, the quantity 
R m i n can be approximated by the following expression, with accuracy to within 1% pQ: 

R min = 0.0488 q~ 0Am A. (23) 

If -R m i n is larger than the magnetosphere radius r m [see (I2"2"j) ]. the magnetic field does not 
exert an important influence on the motion of matter. The flow can circle around the star 
and eventually intersect itself at some point. The further evolution of this flow leads to the 
formation of an accretion disk in the system. If the minimum distance -R m i n is smaller than 
the magnetosphere radius r m , the flow ends up in the zone where the magnetic field exerts 
a substantial influence at some point in its trajectory. The action of the electro-magnetic 
forces in this zone leads to deceleration of the flow and a loss of its angular momentum. As 
a result, the flow is not able to circle around the star and form an accretion disk. Thus, 
the boundary between intermediate polars (which form accretion disks) and polars (which do 
not) is determined by the condition r m = -R mm . Substituting the parameters for the SS Cyg 
system into this relation, we find that the magnetic field separating these two regimes should 
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Figure 10: Same as Fig. |8]for model 6 (B a = 5 x 10 7 G). 

be B a w 10 6 G. However, this estimate of the magnetic field separating intermediate polars 
and polars has a fairly general character, since the values of q vary only slightly in cataclysmic 
variables, and are roughly equal to 0.5. 

Figure [T2] shows the dependence of the accretion rate M a on the magnetic field B a at 
the accretor surface obtained in our computations. The vertical line segments illustrate the 
character of the variations in the accretion rate. The main property of this dependence is its 
non-monotonic character. The accretion rate grows for fields B a < 10 6 G, and the amplitude 
of the accretion-rate variations decrease. A maximum accretion rate is reached at B a = 10 6 G. 
With further increase in B a , the accretion rate falls off. 

This dependence fully corresponds to the reasoning presented above. When B a < 10 6 G, an 
accretion disk forms in the system, and the accretion rate is determined by angular momentum 
transfer in the disk. Increasing the magnetic field leads to an increase in the efficiency of 
magnetic braking in the disk. Therefore, the M a (B a ) dependence acquires a positive slope 
when B a < 10 6 G. Beginning with a magnetic field of B a = 10 6 G, no accretion disk forms in 
the system, and the flow acquires the character of accretion-column streams. In this case, the 
accretion rate M a is determined by the penetrating ability of the stream. The stream cross 
section decreases and the accretion rate falls off with growth in the magnetic field. 
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Figure 11: Same as Fig. |8]for model 7 (B a = 10 8 G). 



4 Conclusion 



We have presented a numerical model designed to investigate mass transfer in semi-detached 
binary systems with strong accretor magnetic fields. Our computations assume the accretor 
magnetic field to be dipolar, with the magnetic axis inclined to the rotational axis. At the base 
of the model is the assumption that the plasma dynamics are determined by the slow mean flow, 
which forms a backdrop for MHD waves that propagate with high velocities. The equations 
describing the slow motion of the matter were obtained by averaging over rapidly propagating 
pulsations. The strong external magnetic field plays the role of an effective fluid with which the 
plasma interacts. In the equation of motion, the mean electro-magnetic force exerted by the 
accretor magnetic field has a form analogous to the friction force between different components 
in a plasma consisting of several types of particle. Magnetic diffusion due to current dissipation 
in turbulent vortices and magnetic buoyancy is included in the numerical model. Moreover, 
averaging the induction equation over rapidly propagating MHD waves leads to additional 
dissipation of the magnetic field (wave diffusion). 

We have used the model we have developed together with a modified, three-dimensional, 
parallel numerical code to simulate MHD flow structures in close-binary systems whose pa- 
rameters correspond to those of SS Cyg. We have presented computational results for various 
magnetic fields B a at the accretor surface. An accretion disk forms in the system in models 
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Figure 12: Dependence of the accretion rate on the magnetic field at the accretor surface B a . 
The vertical line segments show the characteristic amplitudes of the variations of the accretion 
rate. 

with relatively weak magnetic fields B a < 10 6 G. As B a increases, the size of the accretion 
disk decreases and the magnetosphere radius increases. The accretion disk degenerates into 
a spiral disk when B a = 10 6 G, in which the matter is able to undergo only one to two or- 
bital revolutions before falling onto the star. These types of flow apparently correspond to 
the conditions in intermediate polars. No accretion disk forms in models with strong magnetic 
fields B a > 10 6 G. The flow takes the form of accretion-column streams originating from the 
inner Lagrange point L\ and ending at the magnetic poles of the accretor. The morphological 
structure of such flows correspond to polars. 

The field B a = 10 6 G separating these two flow regimes is determined by the relationship 
between the minimum distance to which the stream of matter approaches the accretor i? m i n 
and the magnetosphere radius r m . An accretion disk forms when -R m ; n > r m . Otherwise, the 
flow ends up at some point in the zone near the accretor where the magnetic field exerts a 
substantial influence, preventing the accretion disk from forming. Note that this estimate of 
the magnetic field separating intermediate polars and polars has a fairly general character, since 
it depends only weakly on the parameters of the system. 

Our computational results show that the dependence of the accretion rate on B a is non- 
monotonic. The accretion rate grows for B a < 10 6 G, reaching its maximum value when 
B a = 10 6 G, with the accretion rate falling off upon further increase in B a . The increase in 
the accretion rate for B a < 10 6 G is due to magnetic braking in the accretion disk. When 
B a > 10 6 G, no accretion disk forms in the system, and the flow acquires the character of 
accretion-column streams. Therefore, the stream cross section decreases with growth in B a , 
and the accretion rate falls off in this case. 
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Application 

A Plasma dynamics in a strong magnetic field 

Plasma dynamics in a strong external magnetic field are characterized by the relatively slow 
mean motion of particles along the magnetic lines, their drift across the field lines, and the 
propagation of Alfven and magneto-acoustic waves with very high speeds against the backdrop 
of these slow motions. Over the characteristic dynamical time scale, the MHD waves are able to 
traverse the flow region (along a column-like stream, for example) many times. Therefore, we 
can investigate the mean flow pattern, considering the influence of fast pulsations by analogy 
with wave MHD turbulence. To describe the slow motion of the plasma itself, it is necessary to 
separate out rapidly propagating fluctuations and apply a well defined procedure for averaging 
over the ensemble of wave pulsations. 
Let us consider the relation 

E + -(vxB) = ^, (24) 
c a 

which expresses Ohm's law for the plasma in a magneto-gas-dynamical approximation [27J. 
Here, E is the electric field in the plasma, B = B* + b the total magnetic field, j the current 
density, and a the conductivity. We will represent all dynamical quantities as sums of mean 
values and fluctuations; for example, b = (b) + 5b. Averaging f[2~41) . we find 

(E) + - ((v) x (b» + - ((v) x B*) + - (8v x 5b) = (25) 
c c c a 

The last term on the left-hand side can be estimated using the following expression, often 
applied in dynamo theory (see, for example [3"6"| I3"7]): 

(5vx5b) = a(b)-r] w (Vx(b)), (26) 

where a is the mean helicity of the flow and rj w is the diffusion coefficient for the mean mag- 
netic field due to wave MHD turbulence. We neglect the first term (the a effect), since it 
describes the relatively weak and slow generation of the mean magnetic field in the accretion 
disk. Substituting the resulting expression into (125]) yields 

c(E) + (v) x (b) + (v) x B, - Vw (V x (b» = r/oc (V x (b» , (27) 

where z/od = c 2 / (Aira) is the Ohmic-diffusion coefficient for the magnetic field. 

Let us compare the relative contributions of terms in the resulting expression. Since, as 
a rule, rj w 3> //od, the right-hand side can be neglected. Further, we can neglect the second 
term compared to the third in the case of strong external magnetic fields. Finally, averaging 
Maxwell's electro-magnetic induction equation, we obtain 

Vx<E> = -If>. ( 28 ) 
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Thus, the first term in (1271) is associated with variations in the mean magnetic field (b) over 
the characteristic dynamical time scale. Therefore, this term corresponds to the second term 
to order of magnitude. Substituting the dominating terms into ( 127)) leads the expression 

(v) x B* = r) w (V x (b)) . (29) 

This relation can be used to calculate the mean electro-magnetic forces in the equation of 
motion. Neglecting fluctuations of the density, the wave magnetic pressure and the magnetic 
tension, we find 

(B x (V x B)) = (b) x (V x (b)) + B.x(Vx(b)) ^ 
47rp Airp Anp 

The first term on the right-hand side describes the electro-magnetic force due to the intrinsic 
magnetic field of the plasma, (b) (with the opposite sign). The second term can be calculated 
using (129)) . We then have 

B, x (V x (b)) = B, x ((v) x B») = (yU (31) 
4np 4:7ipr] w r ' 

where (v)j_ is the component of the mean plasma velocity perpendicular to the magnetic field 
B*, and the characteristic relaxation time is 

4.npr} w 



Bl 

B Some estimates 



(32) 



The diffusion coefficient for the magnetic field r] w associated with wave MHD turbulence can 
be estimated using the expression 

^ = y(^ 2 ), (33) 
where r w is the correlation time for pulsations. This expression can be parametrized as follows: 

Vw = ol w — , (34) 
V47rp 

where a w is a dimensionless parameter that is close to unity and l w is the characteristic spatial 
scale for the pulsations, for which we can use the scale for inhomogeneity in the external 
magnetic field, l w = B*/|V.B*|. Since the accretor magnetic field is dipolar in our model, 
L ~ r/3. 

Let us now estimate the dependences of r] w and r on the radius r in a column accretion 
flow. Assuming a uniform density distribution in the flow, we obtain from (1321 . (j34|) rj w oc r -2 , 
r « r 4 . However, the perpendicular cross section of the flow should change with approach 
toward the accretor. We can take this effect into account using the conservation of mass, 
pvS = const. The area of the perpendicular cross section S of the flow can be estimated 
by applying the conservation of magnetic flux, B*S = const. We obtain for a constant flow 
speed v r] w oc r~ 3 / 2 , r oc r 5 / 2 . If we estimate the flow speed as being roughly equal to the 
free-fall velocity, v oc r -1 / 2 , then r/ w oc r -1 , r oc r 3 . Thus, in all cases, the diffusion coefficient 
7] w increases with approach toward the accretor, while the relaxation time r decreases. This 
means that wave MHD turbulence will lead to the virtually complete dissipation of the mean 
magnetic field of the plasma (b) near the accretor, and that the motion of matter will occur 
primarily along the magnetic lines. 



18 



References 

[1] B. Warner, Cataclysmic Variable Stars (Cambridge Univ., Cambridge, 1995). 

[2] A.J. Norton, J.A. Wynn, R.V. Somerscales, Astrophys. J. 614, 349 (2004). 

[3] C.G. Campbell, Magnetohydrodynamics in binary stars (Kluwer Acad., Dordrecht, 1997). 

[4] V.M. Lipunov, Astrophysics of Neutron Stars (Nauka, Moscow, 1987; Springer, Heidelberg, 1992). 

[5] A.G. Zhilkin, D.V. Bisikalo, Astron. Zh. 86, 475 (2009) [Astron. Rep. 53, 436 (2009)]. 

[6] A.G. Zhilkin, D.V. Bisikalo, Advances in Space Research 45, 437-444 (2010). 

[7] A.G. Zhilkin, D.V. Bisikalo, Astron. Zh. 87, 913 (2010) [Astron. Rep. 54, 840 (2010)]. 

[8] T. Tanaka, J. Comp. Phys. Ill, 381 (1994). 

[9] K.G. Powell, P.L. Roe, T.J. Linde, T.I. Gombosi, D.L. De Zeeuw, J. Comp. Phys. 154, 284 
(1999). 

[10] A.G. Kulikovskii, N.V. Pogorelov, and A. Yu. Semenov, Mathematical Aspects of Numerical 
Solutions to Hyperbolic Systems of Equations (Fizmatlit, Moscow, 2001) [in Russian]. 

[11] D.P. Cox, E. Daltabuit, Astrophys. J. 167, 113 (1971). 

[12] A. Dalgarno, R.A. McCray, ARA&A, 375 (1972). 

[13] J.C. Raymond, D.P. Cox, B.W. Smith, Astrophys. J. 204, 290 (1976). 

[14] L. Spitzer, Physical Processes in the Interstellar Medium (Wiley, New York, 1978; Mir, Moscow, 
1981). 

[15] D.V. Bisikalo, A. A. Boyarchuk, P.V. Kaigorodov, O.A. Kuznetsov, Astron. Zh. 80, 879 (2003) 
[Astron. Rep. 47, 809 (2003)]. 

[16] A.G. Zhilkin, Matem. Modelir. 22, 110 (2010). 

[17] D.A. Frank-Kamenetskii, Lectures on Plasma Physics (Atomizdat, Moscow, 1968) [in Russian]. 

[18] F. Chen, Introduction to Plasma Physics (Springer, New York, 1995; Mir, Moscow, 1987). 

[19] B.A. Trubnikov, Plasma Theory (Energoatomizdat, Moscow, 1996) [in Russian]. 

[20] A.R. King, Monthly Notices Roy. Astron. Soc. 261, 144 (1993). 

[21] G.A. Wynn, A.R. King, Monthly Notices Roy. Astron. Soc. 275, 9 (1995). 

[22] G.A. Wynn, A.R. King, K. Home, Monthly Notices Roy. Astron. Soc. 286, 436 (1997). 

[23] A.R. King, G.A. Wynn, Monthly Notices Roy. Astron. Soc. 310, 203 (1999). 

[24] N.R. Ikhsanov, V.V. Neustroev, N.G. Beskrovnaya, Astron. and Astrophys. 421, 1131 (2004). 

[25] A.J. Norton, O.W. Butters, T.L. Parker, G.A. Wynn, Astrophys. J. 672, 524 (2008). 

[26] S.D. Drell, H.M. Foley, M.A. Ruderman, J. Geophys. Res., 70, 3131 (1965). 

[27] L.D. Landau and E.M. Livshitz, Electrodynamics of Continuous Media (Nauka, Moscow, 1982; 
Pergamon, New York, 1984). 

19 



[28] A.G. Zhilkin, Zh. Vychisl. Mat. Mat. Fiz. 47, 1898 (2007) [Comput. Math. Math. Phys. 47, 1819 
(2007)]. 

[29] P.J. Dellar, J. Comp. Phys., 172, 392 (2001). 

[30] A. A. Samarskii, The Theory of Differential Schemes (Nauka, Moscow, 1989; Marcel Dekker, New 
York, 2001). 

[31] A.G. Zhilkin and D.V. Bisikalo, in Proc. 10th Zababakhin Sci. Readings, Intern. Sci. Conf., 
Snezhinsk, Chelyab. Reg., Russia, 15-19 March 2010, www.vniitf.ru/images/zst/2010/sec6/6- 
7.pdf. 

[32] F. Giovannelli, S. Gaudenzi, C. Rossi, A. Piccioni, Acta Astronomica 33, 319 (1983). 

[33] A. A. Boyarchuk, D.V. Bisikalo, O.A. Kuznetsov, VM. Chechetkin, Mass transfer in close binary 
stars (London: Taylor & Francis 2002). 

[34] D.V. Bisikalo and A.M. Fridman, Usp. Fiz. Nauk. 178, 577 (2008) [Phys. Usp. 51, 551 (2008)]. 

[35] S.H. Lubow, F.H. Shu, Astrophys. J. 198, 383 (1975). 

[36] E.N. Parker, Cosmical Magnetic Fields (Claredon Press, Oxford, 1979), p. 608. 

[37] A. A. Ruzmaikin, D.D. Sokoloff, and A.M. Shukurov, Magnetic Fields of Galaxies (Nauka, 
Moscow, 1988; Kluwer, Dordrecht, 1988). 



20 



